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Abstract: In this paper we present a modification of the usual Proteus mirabilis Swarm model. For 
the obtained model (which is a two phase model with a non-linear diffusion term containing memory) we 
set up a collection of a priori estimates. Those estimates allow to get an existence and uniqueness result. 

1 Introduction and results 

Proteus mirabilis is a bacterium that can be either a short cell we call "swimmer" or an elongated 
cell capable of translocation we call "swarmer" . A model of behaviour of Proteus mirabilis colony 
has been proposed by Esipov and Shapiro [5] , based on ideas of Gurtin [TO] . 

In this paper, we prove an existence result to a model which is, in a way, a generalization but also 
a regularization of the Esipov and Shapiro [8] model. 

The model under consideration here is a two phase model with a non-linear diffusion term 
containing memory for one of the two phases. It involves two functions p and Q. The function 
p = p(t,a,x) refers, at time t G [0, T), < T < +oo, to the density of swarmers of age a e [0, A), 
< A < +oo at position x € fi, where 51 is a regular sub-domain of K 2 , with boundary dVL. In 
each point x of dVL, v = v (x) stands for the unit normal vector pointing outside O. The function 
Q = Qipi x ) stands, at time t, for the biomass density of swimmers on f2. 
For a constant r, those two functions are supposed to satisfy the following system: 

^= l —A Q + J^ p(;a,-)e a / T p(;a,-)da + X (A)p(;A,-)e A / T , on [0, T) x Q, (1.1) 

^ + ^ = -HP + V-[{D(M,Q,P) + d)Vp\, on[0,T)x[0,A)xr>, (1.2) 

p(-,0,-) = -Q, on[0,T)xft, (1.3) 

T 

p(p,;-)=Po, on[0,A)xfl, (1.4) 
<2(0,-)=<2o, onfi, (1.5) 

^r=0, on [0, T) x [0, A) x dfl. (1.6) 
o v 

Above, x(A) is an artifice allowing to take into account a possible maximum age A beyond which 
swarmers cannot exist. It has the following definition: 

X (A) = 1, Vie K, x(+oo) = 0. (1.7) 



*This work has been partially supported by "La Region Bretagne", F-35031 Rennes, Program 1042: "Renouvcllc- 
ment des competences dans les laboratoires de recherche", Operation A1C872. 

tLMAM et Lemel, Universite de Bretagne Sud, Centre Yves Coppens, Campus de Tohannic, F-56000, Vannes. 
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Denoting by the space of functions having continuous and bounded derivatives up to order k, 
p = fi{t, a, x) is a function such that 

p 6 Cf([0,T) x [0, A) x fi), /it > 0, lim a, a;) = ~p uniformly in x and t, (1.8) 

with 71 > c(l — x(^))i f° r a constant c > 0. The function £ = £(£, Q) satisfies 

£gC 2 ([0,T) xl), < £ < 1. (1.9) 

In the second equation, V stands for the gradient with respect to the x— variable, and V- for the 
divergence. The diffusion coefficient is the sum of a constant, a priori small, 

d > 0, and of D = D(A4, P, Q) 1 a non negative C\ function of its arguments. (1-10) 

In (|TTT0| . P is defined for < a min < A by 

r A 

P(t,x) = p(t,a,x)e a/T da, (1.11) 

and Q is given by the first equation of the system. 

The memory (or hysteresis) term M. = M\p]{t, x) keeps information on the value of P in the past. 
For four thresholds, P min < p min < p max < P max , with P mm close to p min and P max close to p max , 
Ai is defined as the solution to: 

B -w - p~^ — f p P ~ v T ) H ^ ~ V- ( Pmm ~p P ) H ^ 

Oh ±max Pmax \*max Pmax J Pmin *min \Pmin *min / 

M(0,-)=M o , (1.12) 
with, denoting Pq = P(0, •), 

M a eCl{tt), 0<M <1, M = Q where P < P min and M = 1 where P > P maa; , (1.13) 
and with 

i?r(p) = if p < 0, H r (p) =p if <p < 1 and F r (p) = 1 ifp > 1. (1.14) 

We now turn to the statement of the main result of this paper. 

Theorem 1.1 Under assumptions fO|), fO|) . fOZl)) and f03J), i/p > G L 1 flW' 2 ' 2 nW' 1,4 ([0, A)x 
Q, e a / T dadx) satisfies 



{p ) 2 e 2a/T dx < b / { Po ) 2 e a/T dx, Va G [0,4), 

/ / \V P o\ 4 e 4a ^dadx <b [ [ \V Po \ i e a / T dadx, (1.15) 

/ [ \V Po \ 2 e 2a/T dx <b [ [ \Vp \ 2 e a/r dx, 
Jo Jn Jo Jn 

for a constant b, and if Qq > G L 1 PI VF 2,2 n W^ 1,4 (f2); i/ien, there exists a unique solution 
(Q,p) G L oo (0,T;({L 1 nW 1 > 2 (n))x(L 1 r)W 1 ' 2 ([0,A)xn, e a / T dadx))) to system [T Tp -lT fy coupled 
with hi. 11}) and hi. 12}) . Moreover, Q > and p>0. 
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The precise definitions of the spaces at work in the Theorem are given in the beginning of section [5] 

We now give references where modelling and mathematical methods are developed on age- 
structured population problem: Gurtin and Mac Camy [IT], Marcati [55], Andreasen [U HI 13]; 
possibly with diffusion: Gurtin [TO], Di Blasio and Lamberti [7J, Di Blasio [6], Mac Camy [21], 
Gurtin and Mac Camy [12] , Busenberg and Iannelli [5] , Langlais [13 [18] , Kubo and Langlais [15] , 
Huang |13j and Esipov and Shapiro [8]. For simulation methods we refer to Lopez and Trigiante 
[20] . Milner [24] Kim [14] . Esipov and Shapiro [8], Medvedev, Kapper and Koppel [23], Ayati and 
Dupont [3]. Concerning the biological description of Proteus mirabilis colony behaviour, we refer for 
instance to Rauprich et al [25] . Gue, Dupont, Dufour and Sire [9] and theirs references. 

The paper is organized as follows: in section [2] we present the way to go from the Esipov and 
Shapiro model to system (|l.l|) - (|1.6|l . Then section [3] is devoted to a priori estimates for the solution 
to (|l.l|) - (|1.6p . By a usual procedure consisting in linearizing and passing to limit, we prove the 
Theorem in section [4j 

Acknowledgements: I would like to thank O. Sire for having introduced to me the swarm model 
of Esipov and Shapiro 8] and for stimulating discussions. 

1 would also like to thank F. Granger who, despite stopping its PhD thesis for personal reasons, 
made the first steps towards the result. 

2 Model 

Proteus mirabilis is a pathogenic bacterium of urinary tract that when standing in liquid medium, 
consists in a usual short "swimmer cell" or "swimmer" . When placed on agar medium, if the bac- 
terium density is large enough, it may undergo a differentiation process producing an elongated cell 
with several nuclei called "swarmer cell" or "swarmer" . Those swarmers are capable of translocation 
allowing the bacterial colony to colonise the medium. 

The macroscopic model built by Esipov and Shapiro [8 describes this swarm phenomenon at 
the colony scale. We shall explain this model now. The swarmer behaviour depends on their own 
age. This dependence is taken into account by introducing the age dependent density of swarmers 
p(t, a, x). The link between this density p and the biomass density is a consequence of the fact that 
the mass of each cell is in direct proportion with its length and that the length growth of a swarmer 
is also in direct proportion with its length. Then the biomass density, at time t, of swarmers of age 
a at position x is p(t, a, x)e a l T , where r is the growth rate of the biomass. The first age-depending 
behaviour of the swarmers is that they actively participate in group migration only after an age a m i n . 
Then the definition of the biomass density P of swarmers capable of active translocation is given 
by (|1.11|) . The second age-depending behaviour is that the swarmers dedifferentiate themselves and 
give swimmers. On this topic, Esipov and Shapiro 8 consider two situations. In the first one (Model 
A) , the swarmers dedifferentiate themselves at a given age a max . The second situation consider that 
swarmers may dedifferentiate at each time with a probability l/o (Model B). 

Now, we are able to write the swimmer evolution equation for the biomass density Q. Its evolution 
results from the classical cellular division, with a characteristic time which is the biomass growth rate 
r, subtracting the proportion of bacteria undergoing differentiation and adding the dedifferentiation 
product. In the case of Model A, the evolution equation for Q is then 

^ = —Q + p(-,a max ,-)e a ™*/ T , Q(0,-) = Oo. (2.1) 
at t 

Here, Qq stands for the initial swimmer density and £/r is the fraction of swimmer population to 



3 



produce swarmers. In the case of Model B, the evolution equation is 

dO 1 - £ f* e a / T 

= -Q+ pi; a,-)— da, Q(0,-) = Q o . (2.2) 

dt t J a 

We turn to the evolution of the swarmer density p. Its evolution is linked with ageing and the 
dedifferentiation process, but also to swarm. This last phenomenon is modelled by a non linear 
diffusion term with memory. The evolution equation for p is then 

% + ^ = V'[D(M,Q,P)Vp], (2.3) 
in the case of Model A; and, in the case of Model B, it is 

| + | = -ip + V.[ W Q,P)Vp]. (2.4) 
Both of those equations are equiped with the following initial and boundary conditions: 

(a): p(-,0,') = -Q, (b): p(0, .,•) = (), and (c): %L =0. (2.5) 

r av | an 

The first of those three conditions means the fraction £/r of swimmers undergoing the differentiation 
process produces swarmers of age 0. The initial condition on p means that, at the beginning of the 
process, there is no swarmer. The boundary condition means that swarmer cannot leave the domain 
f2. In the diffusion term V- \D{M., Q, P)Vp\ appearing in (|2.3j) and (|2.4p . and modelling the swarm, 
the diffusion factor D(A4, Q, P) depends on the present value of Q and P but also on the history of 
P. 

The term M then keeps in memory informations concerning the history of the swarmer density. 
Esipov and Shapiro [8] defines A4 as being set to 1 in a given point x if the value of P in x reaches 
a threshold P ma x- Then it remains at the value 1 until the value of P in x reaches another value 
Pram < Pmax- Then they suggest to take 

D(M,Q, P) = D~ M li^A exp (-^-) , (2.6) 

^ ' max ' ^ sat ' 

for given values of D and Q sa t and with, H being such that H(p) = if p < and H(p) = 1 if 
P>0, 

P P P 

lip) = (P~ ~^-)Hip - -_^II1) or 7(p) = p - — — or jip) = p 2 or j(p) = 1. (2.7) 

Pnax ±max ±max 

We also mention that Medvedev, Kapper and Koppel [13] studied this model taking, for a given 
value of k, 

D{M,Q, P) - DiQ, P) = (2.8) 

Now, we explain in what sense the model (|1 . ![) - (| 1 .6|) is a generalization and a regularization of 
the Esipov and Shapiro [8] model. 

First, because of the initial condition (|2.5l b). it is an easy game to see that, at least formally, the 
solution p to (|2.3[) or (|2.4[) satisfies 

pit, a, x) = when a > t. (2.9) 
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Hence we can replace (|2.2I) by 

gn i - c r+°° e a / T 

-£ = q -Q+ / p(-,a,-)—da, Q(O,0 = Q o . (2.10) 

dt t J Q a 

Making this allows to take under consideration, with no loss of consistency, initial data p(0, •, •) that 
are not coming to the more general initial and boundary conditions (|1.3|) - (|1.6|) . 
Secondly, equations (|2.ip and (|2.10|) are particular cases of the general equation (jl.ip with assump- 
tion (|1.8p . The case of equation (|2.ip is recovered setting p, = and A = a max and the case of (|2.10l) 
is recovered setting A = +oo and p = l/a. We also see that (jl.2l) with <i = is a general framework 
inside which (|2.3p and (|2.4p may enter directly. 

Concerning the fraction £ of swimmers to produce swarmers, it seems to depend on experimental 
conditions and to be when the swimmer density is high. It seems then to be reasonable to set that 
£ satisfies (fQ|) . 

The first regularization effect we consider in our model (jl.ip - (jl.6l) consists in adding d > in (|1.2p . 
This may be justified by experimental arguments saying that swarmers always experience a small 
random motion even before and after swarming. 

The second regularization effect which constitutes the most visible modification of the model con- 
cerns the memory term Ai. In order to explain the way to go from the definition of M. by Esipov 
and Shapiro [8] to the definition (| 1 . 12[) . we first notice that the term Ai of Esipov and Shapiro could 
be defined formally, in any x, as the solution to the following equation: 

dM 

-gf(-' x ) = i 1 - M (-> x )) 6 {t/P(t,x)=p ma;c } - M(-,x)6 {t/P{ttX)=Pmm} , M(0,x) = 0, (2.11) 

where &{t/p(t,x)=p ma , x } stands for the Dirac measure in instant t where P(t,x) = P m ax- Of course 
this equation has no real mathematical meaning. But formally, if Ai — and if the value P ma x 
is reached at a given time t, the solution of (|2.1ip experiences a jump +1. When P = P ma x, and 
M. = 1, nothing happens and this is what is needed. In the same way, if at a given time t, P — P m im 
nothing happens if M. = and M. experiences a jump —1 if M. — 1. Now, it is easy to see that the 
right hang side of (|1.12p is nothing but a regularization of the right hand value of (|2.1ip . This is the 
reason why we make this choice to define M. . 

3 A priori estimates 

The key point to get the existence result is a collection of a priori estimates satisfied by (Q, p). They 
are mathematical translations of biological properties. In order to set those estimates comfortably, 
we assume that the assumptions and the conclusions of Theorem 11.11 are satisfied. 

3.1 L 1 and L 2 estimates 

For p > 0, we denote 

\\p(t)\\p= Jjp(t,a,x)\Pe a / T dxda}j '\ \\Q(t)\\ p = (J\Q(t,x)^ dx} '\ (3.1) 

and we have the following estimates saying that the total biomass grows exponentially with a growth 
rate r. 

Lemma 3.1 If the assumptions of Theorem \1.1\ are valid, and if the solution (Q,p) given by this 
same Theorem exists, then it satisfies 

llpWlli + IIQWIIi = (llpolli + IIQolli) e t/r . (3.2) 
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Proof. Multiplying equation (|1.2p by e a / T , integrating then with respect to x and a yields: 

^Ml+lim ( [ pe a / T dx]~ [ p(-,0,-)d!B=-||p||i- I f fJ-pe a/T dxda, (3.3) 

dt a^A \J n J J SI T J J n 

which also reads 
d\\ P \\ 



-+x{A)e A ' T f p(-,A,-)dx- f ^Qdx = -\\p\\i- [ f ppe a / T dxda. 
Jn Jn T T Jo Jn 



dt 

Integrating now with respect to x, we obtain 



(3.4) 



o J si 



dt Jn t 
Summing up (|3.4[) and (|3 . 5|) finally gives 

4|M|i + ||Q||i) 1 



ppe a/T dxda + X (A)e A/T [ p(-,A,-)dx. (3.5) 
Jn 



dt t 

proving the Lemma. 



IpIIi + IIQIIi), (3-6) 



The second Lemma concerns the L 2 — norms of p and Q. It mathematically translates that the 
biomass cannot be so gathered that a null area set contains a positive biomass quantity. 

Lemma 3.2 If the assumptions of Theorem \1.1\ are realized and if the solution (Q,p) exists, then 
for any t G [0,T), it satisfies 



|2+||QW||2<c(||p || 2 +||Qo||2), (3.7) 

for a constant c (c only depends on A, t, T and p, — s\xp{p{t 1 a, x), t G [0, T), a G [0, A), x G 0,}.) 
Proof. First, integrating (|1.2j) with respect to x, we get: 



d_ d_ 

dt da 

Making the same, after multiplying (|1.2[) by p, gives: 

l{ll + l)Un ^ dX ) = - In ^ ^ -Jn iD + ^ ** " ° 
Defining, for a fixed a and for p = 1 or 2, 

pP{s,a + s,x) dx, (3.10) 



pdx\ = ppdx < 0. (3-8) 



(3.9) 



we have 



-«•>- ((i + 5)(/ '*))<*'' ,+ "* a 



(3.11) 



Now, since for fixed t and a and with a — a — t we have J n p p (t, a, x) dx = r p (t), if i < a the relation 
r P (t) < »>(0) = J n p p (0, a — t,x) dx reads 

p(t,a,x)dx< / po(a — t,x)dx, / p 2 (t,a,x)dx< / p 2 (a — t,x)dx. (3.12) 

Jsi Jn 
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In the case when t > a, the relation r p (t) < r p (t — a) — J n p p (t — a, 0, x) dx gives 



p(t,a,x)dx< [ —Q(t — a, x) dxda, j p 2 (t,a,x)dx< ( ^jQ 2 (t — a, x) dxda, (3.13) 
in Jn T Jn Jn T 

Secondly, we multiply equation (|1.2p by e a l r p and we integrate in x and a. Since 



da' 



^pe a ' T da = -l ^ p e a ' T da-- I p 2 e a ' T da + 



da 



p 2 e a ^ 



a=0 



and since 



we obtain 



44 

dt 



x(A)e 



A/t 



f p 2 (t,0,x)dx = f ^Q 2 dx, 
Jn Jn T 

[ p 2 (-,A,-)dx + 2 [ [ (D + d)\Vp\ 2 e a/T dxda 
Jn J o Jn 



(3.14) 
(3.15) 



2 f f pp 2 e a ' T dxda = f L I Q 2 dx + -\ 
Jo Jn Jn T T 



(3.16) 

As the second, third and fourth terms of the left hand side of equation (|3 . 16[) are non negative, we 
may deduce 

^<-^IIQII2 + -IIpII2. (3-17) 

dt T Z T 

In the same way, multiplying now (jl.ip by Q and integrating yields 

r A 

-Q 2 dx + 2 



dt 



ppQ e a/T dxda + 2 X {A)e A/r / p(-,A, -)Q dx. (3.18) 
o Jn Jn 



Concerning the second term of the right hand side of the last equality, using Young's inequality and 
formula (|3.12p and (|3.13[) , we get 

/ f ppQ e a/T dxda < [ ( [ p 2 p 2 dx] ( [ Q 2 dx] e a/r da 
Jo Jn Jo \Jn ) \Jn J 

J (J ^Q 2 (t-a 1 -)dxy /2 e a/T da 



oT/t 



< /imax(e /r , — — ) 



+ H(A - t)e l / T J ( J p 2 {a - t, -) e 2(a "*'/ T dx] 
( / Q 2 (t~a, -)dx] 1,2 da 



1/2 



da 



o x Jn 



WQh (3.19) 



+ H(A-t) f ( f pg(a - t, .) e 2 ( Q -*)/ T dx] 



1/2 



da 



WQh, 



where H(a) — if a < and H(a) = 1 if a € [0, +oo] and where p = sup{p(t, a,x), t 6 [0, T), a <E 

[o,A),xen}. 

The third term of the right hand side of (|3.18p may also be estimated: 



X (A)e A/T / p(;A,-)Qdx< X (A)e A ^[ / p 2 (;A,-)dx 
Jn \Jn 



1/2 



Q 2 dx 



1/2 



(3.20) 



Hence, applying again (|3.12[) and (|3.13p . 

1/2 



x(A)e A /* f p(-,A, -)Q dx < X {A)e A / T (f Q 2 (t-A, •) da) \\Q\\ 2 , if t 

/ r \ 1/2 

< X {A)e 2t ' T {J^p 2 {A-t r )e 2 ^' T dx\ \\Q\\ 2 , if t < A 



Using (|3.19|) . ()3.21j) and ()1.15|) in (|3. 18|) . for four non negative constant c%, c 2 , C3 and C4 we have 

^iinii 2 / /■'//■ \ 1/2 / r \ 1/2 



dt 
Setting 



< fci||Q||2+^ ( y g 2 (t-a,-)^) 1/ Wc3|| i oo||2+c 4 (ff(<-A) y Q 2 (t-A,-)dx 



IQII2. 
(3.22) 



F(t)= sup ||g(s)|||+ sup ||p( S )|| 2 , (3.23) 
se[o,T) se[o,r) 



we have 

F'(t) < max(0, «M ) + m ax(0,^), (3.24) 
and from (j3.17p and (I3.22[) we deduce 

F'(t) < {( Cl + c 2 T + c 3 + c 4 )^/FjF))^F(t) + + ^)F(t) < c 5 F(t), (3.25) 
for a constant C5, from which we deduce that 

F(t) < e CsT F(0), (3.26) 
and, as a consequence, that for a constant C6, 

IIQWIIi+llpWIIl^CBdlQolll + IIPolll). (3-27) 

Finally we get f|3 . T|) as a consequence of (|3.27p , this ends the proof. ■ 

As a consequence of the Lemma 13.21 we have the following Corollary whose biological meaning 
is: Biomass cannot be created from nothing. 



Corollary 3.3 If the assumptions of Theorem \l.l\ are realized and if po = and Qq — then the 
solution (Q,p) given by the Theorem satisfies 

p = and Q = 0. (3.28) 

In order to establish the previous estimates, we have assumed that Q > and p > 0. We can show 
that this is a consequence of the non negativity of Qo and po. 

Lemma 3.4 If the assumptions of Theorem \l.l\ are true and if there exists a solution (Q, p) to system 
iTO P -fO p . then Q > and p > 0. 



The proof of this Lemma is close to the one of Lemma [ 
Proof. We define 

p-=min(p,0), Q~ =min(Q,0). (3.29) 



Of course, 



p-(Q,;-) = 0, Q-(0,-) = 0, p-(-,0,.) = -Q- (3.30) 

T 
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Now, multiplying (|3.30[) by e a / T p and integrating, since 



dp - dp 

m p ~ = -m p ~' (3 ' 31) 



proceeding as while establishing (|3.17p . we get 



d\\p-\\i / 1 M „_ lia , i 



dt - r 2 "" "* ' r 
Multiplying now (jl.ip by Q~ and integrating gives 

d\\Q-\\l _ f l-Z„-2 



<Zs\\Q-\\2 + z\\P~\\i (3-32) 



dt 



= 2 f - — ^Q- 2 dx + 2 f f ppQ-e a/T dxda + 2 x (A)e A/T f p(-,A,-)Q- dx 
Jq T Jo Jq Jq 

rA 



<2 f - — -Q- 2 dx + 2 [ [ pp-Q-e a/T dxda + 2 X (A)e A/T [ p~(-,A, -)Q~ dx. 
Jq t Jo Jq Jq 



(3.33) 



lo Jq 

Since, as we had (|3.12[) and (|3. 13[) . we have here: 

f -2 f _2 f _2 f £ 2 _2 

/ p dx — p dx if t < a and / p dx = ~Q (t — °j ') c^^c if t > a, (3.34) 
we can finish the proof as in the proof of Lemma 13.21 and get 

||Q-(*)|| + < c'(\\Qo\\ + \\Po II) - 0, (3.35) 

giving the Lemma. ■ 

3.2 Estimates on the derivatives 

As a by product of the proof of Lemma 13.21 we can deduce from (|3.16[) the following result insuring 
a first control on the regularity of p. 



COROLLARY 3.5 If the assumptions and the conclusion of Theorem IJ.il are valid then, for any 
< s < T, 



\\Vp\\%dt= [ [ [ \Vp\ldxdadt < 3, (3.36) 
Jo Jo Jq 



I I Vol 2 d rdrull < 

d 

for a constant c (depending only on A, t, T, p, \\p0W2 and \\Q0W2) 
Proof. Integrating (|3.16p from to s yields: 

rA 



x(A)e 



A/t 



f f p 2 (-,A, -)dxdt + 2 fff (D + d)\Vp\ 2 e a ' T dxdadt 
Jo Jq Jo Jo Jq 



2 I I [ PP 2 e a/T dxdadt= [ f ^Q 2 dxdt+- I \\p\\* dt + ||p(0)||l - \\p{s)\\% (3.37) 
Jo Jo Jq Jo Jq t t 



Using the previous estimate concerning p and Q, we may deduce 



d f f f \Vp\ 2 e a/T dxdadt < c, 
Jo Jo Jq 



(3.38) 



and (|3. 36[) follows. 
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Because of the form of the non linearity in (jl.ip ~ p.6p . we need a supplementary estimate con- 
cerning 



/ ||Ajo|||dt= f [ I \Ap\ldxdadt. (3.39) 
Jo ' Jo Jo Jn 

This estimate is a consequence of an estimate on || Vp||4 and on 

||p|| o = sup{|/>(i,a,a:)|,f G [0,T),ae [0,A),x G fi}, ||Q||oo = sup { \Q(t, x)\, t G [0,T),a;Gfi}, 

(3.40) 

that we now set. 

Lemma 3.6 The solution (Q,p) given by Theorem \l.l\ satisfies 

||p||oo + ||Q||oc<fc, (3.41) 

where k is a constant depending only on T , A, sup te [ 0T ) ||p(t)[[2 ^nd sup tg j T ) HQWIh (which are 
estimated by Lemma\37~ 



Proof. As we already see that ||p(i)||2 is bounded, using a method similar to Ladyzenskaja, 
Solonnikov and Ural'ceva [TOj, (paragraph III-8) we deduce that 

\x(A)p(;A-)e A / T \ <k x , (3.42) 

where the constant k\ only depends only on A and sup te m r ) |p(t)||2- Then defining 

P = / pe a/T da, (3.43) 



we deduce from (jl.2p that P is solution to the following parabolic equation 
dP 



{A)p(; A, -)e A ' T = - f ppe a ' T da + V-({D + d)VP) + -Q, (3.44) 
Jo T 

from which we get that 

||P|| oo =Hup{|P(t,a;)|,te[0,T),a:en}<A 2) (3.45) 

f A 

where ki only depends on sup ||p(i)||2, sup ||Q(t)||2 and sup |Po| = sup | / poe a ^ T da\ which is 

te[o,T) te[o,r) xen ieo Jo 

finite by assumption. 

Then, (|3.45p and (|3.42p give that ||Q||oo is finite, and as a consequence, (|1.2p and (|1.3p finally give 
the bound on ||/o||oo7 ending the proof. ■ 



Lemma 3.7 The solution (Q,p) given by Theorem \l.l\ satisfies 

||VQ||4+||Vp|| 4 = (jjVQtdx^J ' + JjVp\ 4 e a ^dxda^j 1 <^, (3.46) 

for a constant C (C does not depend on d). 
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Proof. Multiplying equation (|1.2p by —V • (| Vp| 2 V / o)e a / T , and integrating in a and x gives: 

- f I V • (|Vp| 2 V/>) e a/T dxda -I f |^ V • (|V^| 2 Vp) e a/r dxda = 

/ [ ppV ■ (\Vp\ 2 \7p) e a/r dxda- [ ( V ■ [(P + d)V p]V ■ {]V p\ 2 V p) e a/r dxda. (3.47) 
Jo Jo Jn 

Making a double integration by part, and following a straightforward computation procedure, we 
get 



f f V ■[(D + d)Vp]V-(\Vp\ 2 Vp)e a/T dxda = 
Jo Jn 

[ [ (D + d)(\Vp\ 2 \V 2 p\ 2 + 2H(Vp,V 2 p)) e a / T dxda + E, (3.48) 
Jo Jn 



where 



l^ = (0) 2 + ^> 2 + (gf, (3.43) 

H , y V 2 ^f dpd 2 p | dp <9 2 p \ 2 | <9 2 p | dp <9 2 p \ 2 

V9xi 9x 2 dxidx\dxi) \dx\dx\dxi dxidx\) 



and 

-tit 



dP(M,Q,P) dp 



2 d2 P ^ 2 ( dp I dp ) dp )c a/T dxda (3 51) 



Vp| 

dxidxj dxi dxidxj 8x2 dxjdx-i dx 



Since 

dD(M,Q,P) _ dP dM 8P_dQ_ 3D dP 

dxi dM dxi dQ dxi dP dxi ' 



(3.52) 



in view of the regularity of P, of equation (|1.12[) that gives a control on dM jdxi in terms of dPj 9a;,, 
we get for a constant C\ 



|E| < Ci / f (|VP| + VQ|)|Vp| 3 \V 2 p\e a / T dxda 
Jo Jn 



<lj J l V -°| 2 l V V| 2 e a ' T dxda + ^ / / (l VF l + l v QI) 2 | v H 2 e a / T dxda; (3.53) 

in order to get the last expression in (I3~53l . we used f/V < ft/ 2 + ±V 2 with V = Ci(|VP| + 
|VQ|)|Vp|. 

£^ I I 4: ^1^7 1^ 

Concerning the other terms of (|3.47p . since — — — = 2 — — — |Vp| 2 = 4 ^ ■ Vp |Vp| 2 , making an 
integration by parts, we get 

" f f t V-(\Vp\ 2 Vp)e^dxda = \ d -^t. (3.54) 
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In a similar way, 



d{\Vp\ A e a ' T ) , , 1 „_ ... 



\ X {A) J \Vp(; A, -)| 4 e A ^dx --^zj (e 4 + ( J|q) 4 ) I VQ| 4 - ^ |] Vp|| 4 , (3.55) 



and 



C 2 (||Vp|| 4 + l), (3.56) 



f [ pp V • (|Vp| 2 Vp) e a ' T dxda 
Jo Jn 

= - [ [ p\Vp\ 4 e a/T dxda- [ [ V p p\V p\ 2 e a/r dxda 
Jo Jn Jo Jn 

using the regularity of p and the estimate on sup(p) given by Lemma [J 

The regularity of £ and the estimate on sup(Q) give |£ 4 + (§^Q) 4 \ < C*3 for a constant C3. Hence 
(gUIII yields 

I^T^ + 4 x(A) / |Vp(-,A0| 4 e A/T d*+ / A / (£ + d)(|Vp| 2 |V 2 p| 2 ) e^dsda 

+ 2/ / (Z? + d)H(Vp,V 2 p)e^dx<^ / |VQ| 4 dz + i||Vp||t + C< 2 (||Vp|| 4 + l) 



d 

4 



/ / \Vp\ 2 \V 2 p\ 2 e a/T dxda + 
Jo Jn. 



d 



|Vp| 4 e a ' T dxda + C 4 (\\VQ\\i + || Vp|| 4 ) ) , (3.57) 



10 Jn 

and passing the fourth term of the right hand side in the left hand side we can deduce 

— j t — < -^-(l|vp|| 4 + ||VQ|| 4 + l), 



(3.58) 



for a constant C5. 

Multiplying equation (|1.2|) by — V • (|Vp| 2 Vp)e 4a / r , and making the same operations as previ- 
ously, we obtain an inequality which is (|3.57p with e a l T replaced by e 4a / T and ||VQ||| replaced by 
Jo In \ V P\ 4 e 4a/T dxda. From this, we can deduce 



d(I A f a \Vp\*e^dxda) 



(lt 



+ X (A) / \yp(-,A,-)\ 4 e 4A ^dx 



< 



C 6 



[ [ \Vp\ 4 e 4a/T dxda+ ||VQ|| 4 + l). (3.59) 
Jo Jn ' 



On another hand, computing the gradient of (jl.ip . and multiplying by VQ|VQ| 2 yields 



id||vg|| 4 1 r da 



-Q\VQ\ 4 dx 



|VQ| 4 dx 



n \jq 



Wppe a/T da) ■ VQ\VQ\ 2 dx 



4 dt t J n dQ 

+ Is (/ Wpe^dd^ -VQ\VQ\ 2 dx + x (A) £ (v j o(-,A,-)e A/T ) -VQ\VQ\ 2 dx. (3.60) 
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Because of the regularity of £ and p and of Lemma 13.61 since, applying Young's inequality, 



J (J^ pVpe a/T da) -VQ\VQ\ 2 dx 

~ (/ {Jo ^ P ^ a/Tda ) da; ) (/ (l V ^l 3 ) ' dx 



3/4 



<»(j £\Vp\ 4 e 4a / T dxda\ \\WQ\\l, (3.61) 



and 



X (A) [ {Vp(-.A, •) e A ^) . VQIVQI 2 dx < X (A) ( [ |Vp(-, A, -)| 4 e 4A ^ dx) \\VQ\\l 
Jn \Jn / 

< x(A)^f ([[ jT* J | V P | 4 e^dxda) ^ + ||VQ|| 4 + l) ||VQ||t, (3.62) 



we deduce from (|3. 601) 
d||VQ|| 



^<c 4( ||vg|| 4 + i) + f 



^ y |Vp| 4 e 4a / r do;da) 1/4 + ||Vg|| 4 + l^)||VQ||i (3.63) 



Inequalities (|3.58p . (|3.59[) , (|3.63|) and the assumptions on Q and /9 gi ye 

rA 
10 JQ 



1/4 C 



II vqiu + II vp|| 4 + ( jf / o |Vp| 4 e 4a /^d a ) 

and finally the Lemma. 

Lemma 3.8 TTie solution (Q,p) given by Theorem \l.l\ satisfies 



l|VQ(t)||i<^, ||Vp(*)||l<^ foranyO<t<T, 



(3.64) 



f \\Ap(t)f 2 dt= f /" / |Ap(t)| 2 e Q /^a^<-J, /, 
o Jo Jo Jn a 



or any < s < T, 



(3.65) 
(3.66) 



dt [ ) 



dt < c, 



dt {) 



dt < —r anc 
> a- 



dt < — 7t for any < s < T, 



(3.67) 



for a constant c (which does not depend on d). 



13 



Proof. Multiplying (|1.2|) by (— A j o)e a / T and integrating in a and x gives 
dt 



d\\Vp\\l , 



(A)e A / r / |Vp(-,A,-)| 2 rfa; + 2 / / /x| Vp| 2 e a / T dxda 
Jn Jo Jn 

+ 2 [ [ (D + d)\Ap\ 2 e a/T dxda= [ [ V > ■ V p pe a/7 dxda 
Jo Jn Jo Jn 

- f f \7[D(M,Q,P)]-\7pApe a / r dxda+-\\Vp\\l+ f + (§-Q) 2 )\VQ\ 2 dx 
Jo Jn > Jn 

< Cl (||Vp|| 2 + ||V P ||2 + ||VQ||2) 
< Cl (||Vp|| 2 + ||Vp||| + ||VQ||i) + iy" 



o Jn 

A 







clt> 


(l>- 




dD 
+ dQ 










+ 4, 


o Jo 



for a constant ci. Now transferring the last term of the right hand side in the left hand side and 
using the estimate on ||V/o||4 and ||VQ||4 which also give an estimate on ||V.M||4 and ||VP||4, we 
get, for a constant C2 

+ ^\\Ap\\j < c 2 (\\Vp\\t + \\VQg + ~), (3.69) 

In a similar way, we can also get 

d(J A J Q \Vp\ 2 e^dxda) , r A 
dt 

and 



d||VQ|' 2 



2 



dt 



~ C4 Uo Is \ V P\ 2 e2a/Tdxda+ ( 3J1 ) 



From the three last inequalities we get (|3.65p . Integrating (13.69)) from to s gives p.66[) . The 
estimate on ||VQ(t)||2 is then obtained as the estimate on ||VQ(f)||4. Estimate (I3.67P is finally a 
direct consequence of equations (jl.lj) and (|1.2p . ■ 

Remark 3.1 We could also prove that ||V<3||oo and ||Vp||oo are bounded. 



4 Existence and uniqueness of the solution 

Once the a priori estimates are set, the proof of existence is classical and in the spirit of Ladyzenskaja, 
Solonnikov and Ural'ceva [16] . It essentially consists in linearization and passing to the limit using 
the estimates. 

In the following L p (fl) and L p ([0, A) x Q, e a / T dadx) are the functional spaces associated with the 
norms (|3.ip . W k,p (Q) and VF fe,p ([0, A) x f2, e a l T dadx) are the Sobolev spaces composed of functions 
whose derivatives up to order k are in L p (fl) or L p ([0,A) x Q, e a/T dadx). L p ([0,T) x fi) and 
L p ([0,T) x [0,A) x fi, e a / T dadxdt) are the spaces of functions having finite norm 

/ rT pA ? \ i/p [ r T f \ 1/p 

[ill \p(t,a,x)\ p e a/T dxdadt) or / / \Q(t, x)\ p dxdt) , (4.1) 
\Jo Jo Jn ) \Jo Jn J 

and W k - p ([0,T) x Q) and W k ' p ([0,T) x [0, A) x O, e a / T dadxdt) are their associated Sobolev spaces. 
Finally, for a functional space W, L°°(0,T; W) stands for the functions whose norm in W is finite 
for any t G [0,T]. 
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4.1 Linearization 

We linearize the system (jl . 1 [) - (| 1 Then using classical results on pde and ode, we give an existence 

and uniqueness result for the solution to this linearized system. 

We set Q° = Qo and p° = po and for n £ N*, we consider (Q n ,p n ) solution to: 

^1 = Iz1q» + f a, -) e a/ >(-, a, •) da + X (A)p n (; A, -)e A l\ on [0, T) x f2, (4.2) 

r)n n 8n n 

-j^ + -^ = -pp n + V-[(D{M n -\Q n -\P n - 1 ) + d)V P n ], on[0,T)x[0,A)xO, (4.3) 

p n (;0,-) = -Q n , on [0,T)xO, (4.4) 

T 

p"(0, v )=A>, on[0,A)xn, (4.5) 

Q"(0,-) = Qo, onH, (4.6) 

= 0, on [0, T) x [0, A) x 90, (4.7) 

a v 

where for n E N, 



P n (t,x) = / p n (t,a,x)e a / T da, 

** Q>min 



and A4 n is solution to 



dt P-max Prnax \ Pmax Pma 



' g r f Pmm J" ) H r (M n ) 



(4.8) 



Pmin Pmin \ Pmin Pn 

M n (0,-)=M Q . (4.9) 



Theorem 4.1 Under assumptions fOD.fOU. rPZW and /TQ3)) . z/po > £ L 1 nM /2,2 nW /1 ' 4 ([0, A)x 
ft, e a / T dadx) satisfies l[TJ5\) and if Q > £ L 1 n V^ 2 ' 2 n W /1 ' 4 (0) i/ien /or any n £ N, £/iere exists 
a umgue solution (Q",p n ) £ C£(0, T; (L 1 n W x ' 2 (0)) x (L 1 n T^ 1 - 2 ([0, A) x ft, e a / T dadx))) to systom 
( |4 - j^-Tp coupled with |^.&[ ) an d 14 -Sty - Moreover, Q n > 0, p n > and (Q n ,p n ) satisfies estimates 
{kfy . (S7!\) , {gggP , fSSgj, (ESSj and P^7| j mf/i constants independent of n. 

The proof of this Theorem uses only classical pde and ode arguments. Hence we only sketch it. 
Proof. The proof consists in an induction procedure. Because of the assumptions on (Qo, Po) and 
the definition of (Q , p ), the Theorem is true for ri = 0. 

Then, if the Theorem is true for n — 1, by regularization arguments, we can get that that A^ n_1 
exists and is unique on [0, T) x ft and that X™- 1 £ C°(0, T; H/ 1 ' 2 ^)) n C^O, T; L°°(ft)). Hence we 
deduce that, for each I £ N*, there exists a unique solution (Q n ' l ,p n ' 1 ) £ C b °(0,T; (L 1 n W^^ft)) x 
(i 1 nVF 1 ' 2 ([0,A) x ft, e Q / T dada;))) to 

^ = I Q n,i _ |q»,I-i + I a .)e°/>(-, a, •) da + X (^)p n ''(-, A -)e A/T , (4.10) 

dt r t Jo 

9p n,/ Q p n 



dt da 



-pp n ' l + V- [(D(M n ~ 1 ,Q n -\P n - 1 ) + d)Vp n ' 1 ], (4.11) 

p n > l (;Q,-) = ~Q n ' l -\ p n ' l (0,;-)=Po, ^=r, =0, (4.12) 
t a z/ an 

Q"' i (0,-) = <9o ) (4.13) 
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where Q n, ° is defined as Q n '° — Q n 1 . This deduction involves first a classical semi-group or 
Galerkin routine in order to deduce that there exists a unique solution p n < 1 to P~TT|) - gUg) as soon 
as the existence of (Q n,l ~ l , p n,l ~ l ) is achieved. These routines are explained in Lions and Magenes 
[19] . Ladyzenskaja, Solonnikov and Ural'ceva [Hj, or -in a context close to our- in Langlais [I7| . 
Once the existence of p n ' is established, the existence and uniqueness of Q n,t follows. Now, following 
the way leading to (|3~T7f and ([3~22]) we deduce that {Q n ' l ,p n ' 1 ) satisfies 

i_|| Q n,<-l||| +111^11, (4.14) 

d|IQ ^'' lli < (ci ||Q' U || 2 + cillQ^" 1 ^ jf* ( f n (Q n ' l ) 2 (t ~ a, •) dx) V2 da + caHpolla 

+ Ci (H{t - A) f n (Q n ' l ) 2 (t - A ■) dx) 1/2 ^j \\Q n ' l \\ 2 , (4.15) 

which is enough to deduce that (||p n, '||2 + \\Q n H2) is bounded. As a consequence of this bound, we 
get that, for a subsequence still denoted I, {Q n '\p n < 1 ) — ► (Q n ,p n ) in L°° (0, T; (L 2 (ft)) x (L 2 ([0, A) x 
Q, e a / T dadx)) weakly—*, where (Q n ,p n ) is solution to (|4.2[) - (|4.7|) . Finally, the estimates are led 
in the same way as in section [3] The uniqueness follows directly from the linear character of (|4. 2|) - 
(|4.7p . Hence, the Theorem is true for n. 

The induction procedure is then straightforward to end the proof of the Theorem. ■ 



4.2 Existence 

From estimates (|3.65|) . (|3.66|) and (13.67| we deduce that the sequence (Q n ,p n ) is bounded in 
W 1,2 ([0,T) x CI) x W /1 ' 2 ([0, T) x [0, A) x CI). Hence, up to a subsequence still denoted n, we have 
(Q n ,p n ) — > (Q,p) in W^ 2 ([0,T) x Cl) x VF 1 - 2 ([0,T) x [0,A) x O, e a / T dadxdt) weakly, and then, in 
L 2 ([0,T) x O) x L 2 ([0,T) x [0, A) x 17, e a l T dadxdt) strongly. 

From this we can also deduce that P n — > P in L 2 ([0,T) x Cl) strongly, with P defined from p by 
(fTTTTT) . In view of ([4~9f . we can deduce that (X"), (dM"/dt) and, taking the gradient of (|4~9|) . 
(V.M n ) are bounded in L 2 ([0,T) x 17). Then extracting again a subsequence still denoted n, we 
deduce .M™ — > strongly, where .M is solution to p,12[) . 

Using now the regularity of P>, we obtain D(M n - x , <5 n_1 , P"" 1 ) — ► D(M,Q, P) in i 2 ([0, T) x ft) 
strongly. 

The regularity of trace operators gives x(A)p ra (-, A, •) — ► x{A)p(-, A, •), p"(-,0,-) — ► p(-,0, •), 
p n (0, v ) — -> p(0, v ), Q"(0,-) — ► Q"(0,-) weakly, and using (JSH), dp n /dl? lan — > dp/ff^\ga 
weakly. 

Then passing to the limit in (|4.2p - (|4.9[) . we obtain that (Q,p) is solution to (|1.1[) - (|1.6p coupled 
with (fTTTTj) and ([TTL?]) . 

Once this existence result is established, using regularizations and truncations, we can start the 
computations of section [3] giving the additional regularity and the non negativity of Q and p. 
It now remains to prove the uniqueness of the solution. 
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4.3 Uniqueness 

Consider (Q,p) with associated P and Ai and (Q,p) with associated P and J\4 two solutions of 
- (|1.6p . They both satisfy the estimates and the difference (Q,p) = (Q — Q,p — p) satisfies 

S = —Q + f A p(; a, -)e a/ V(-, a, •) da + X (A)p(-, A, -)e A '\ (4.16) 



dt t 

^ + ^ = -pp + V-[(D(M,Q,P) + d)\7p]-V-[{D(M,Q,P)-D(M,Q,P))\7p], (4.17) 

p(;O r ) = -Q, p(0,-,-) = 0, Q(0,-) = 0, |4 =0. (4.18) 
Multiplying (|4.17|) by pe a / r and integrating gives 



+ X (^)e A/T / A ■) dx + 2 / / (D(M,Q, P) + d)|V/5| 2 e Q / T dxda 



n Jo Jo 



= -2 / / MP 2 e a/T dxda + f L-Q 2 dx+-\\p\\ 
Jo Jn Jn T T 

+ f f (D(M,Q,P) - D{M,Q,P))Vp-\7pe a/T dxda 
Jo Jn 

< h(\\p\\i + IIQIIl) + fc 2 ||Vp|| 0O (||p|| 2 + ||Q|| 2 )||Vp|| 2 

< fcidlPlll + IIQIIl) + f \Wp\L(\\ph + WQh) 2 + - a \W~p\\1 (4.i9) 

for constants k\ and k%- Passing the term 4||Vp||| in the left hand side gives, for a constant £3 



d\\~ P \\l / k 



<^{\\P\\1 + WQWD- ( 4 -20) 



dt ~ d 

Making the same, but multiplying (14. 17[) by pe 2a / T yields 
d(J A J n \p\ 2 e 2 ^dxda) /T 



df , x{A)e 2A ' T Jj(;A,-)dx<'^n A Jjp\ 2 e 2a / T dxda+\\Q\\iy (4.21) 



Multiplying (|4.16|) by Q gives 



^M<^±( [ A I irf^Tdnla^WOft). (4.22. 



Jn 



for a constant k^ From the last three inequalities we deduce, for a constant K 



\i + \\p\\t 



[ [ \p\ 2 e 2a ' T dxda< iff||Q| t =o||i + ||p| t =olli+ / / \p lt =o\ 2 e 2a/T dxda) = 0, 
Jo Jn \ Jo Jn / 

(4.23) 

giving Q = p = and then the uniqueness of the solution to (11. 1|) - (|1.6p . ■ 
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